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Abstract 


Monthly maps of sea surface height are constructed for the North Atlantic ocean us- 
ing TOPEX/Poseidon altimeter data. Mean sea surface height is reconstructed using a 
weighted combination of historical, hydrographic data and a synthetic mean obtained by 
fitting a Gaussian model of the Gulf Stream jet to altimeter data. The resultant mean 
shows increased resolution over the hydrographic mean, and incorporates recirculation 
information that is absent in the synthetic mean. Monthly maps, obtained by adding 
the mean field to altimeter sea surface height residuals, are used to derive a set of zonal 
indices that describe the annual cycle of meandering as well ai position and strength of 
the Gulf Stream. 
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2 RECONSTRUCTING THE MEAN SEA SURFACE HEIGHT 


1 Introduction 

Tbit report describe* the construction of monthly maps of total tea surface height {ath) 
In the North Atlantic ocean using TOPEX/Potaidon altimeter data. These maps are 
then used to derive sonal indices that describe the transport, position, velocity, and 
curvature of the Gulf Stream. 

Collinear analysis of the subtrack data is described extensively in [11. The end- 
product of this analysis is cleaned residual ath profiles splined along-track on to a com- 
mon latitude grid. Since the mean is removed along with the geoid in the collinear 
analysis, it is first necessary to recover the mean height field in order to measure total 
aak difference. The first section of this report deals with the reconstruction of a mean 
for the North Atlantic using both historical hydrographic data and a synthetic mean 
derived using the Kelly/ Gille/Qiu method [4], [5]. This mean is used along with the Sah 
residuals to obtain monthly maps of sah. The second section describes the calculation 
of zonal indices that characterize the state of the Gulf Stream (GS) from these monthly 
maps. Contour plots of the monthly maps are included in the Appendix. 

2 Reconstructing the Mean Sea Surface Height 

For the purposes of this study, the aubtracks that fell between latitudes 20°N to 60°N 
and longitudes 30°W to 80°W were analyzed for cycles 4 through 80 which spanned 
a time period from October 1992 to November 1994. Cycles one through three were 
discarded because of poor data quality. Figures 1 and 2 show the coverage of ascending 
and descending ground tracks in this area along with the tah variance along the tracks. 
The peak of the variance reveals the position of the GS jet. 

Two separate calculations of the mean are described in the following subsections. 

2.1 Obtaining the hydrographic mean 

To obtain a climatological mean, we used the Lozier/O wens/ Curry HydroBase data set 
which calculates the mean dynamic height by gridding raw, historical data for pressure, 
temperature and salinity on density surfaces [3]. This yields net CDF files which are 
then used to extract the dynamic height of the surface relative to a reference depth of 
2900 meters. Figure 3 shows a contour plot of these data. 

2.2 Obtaining a synthetic mean 

To generate a synthetic mean using the Kelly /Gille/Qiu method, the Gulf Stream ve- 
locity was modeled as a Gaussian jet *4[. Assuming that the GS meanders a distance 
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2 RECONSTRUCTING THE MEAN SEA SURFACE HEIGHT 
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Figure 1: TOPEX/Postidon ascending ground track coverage for the North Atlantic 
showing tih variance along track. 
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Figure 2: TOPEX/Pottidon descending ground track coverage for the North Atlantic 
showing ash variance along track. 
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2 RECONSTRUCTING THE MEAN SEA SURFACE HEIGHT 
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Figure 3: Dynamic height from historical hydrographic data, reference level 2000 meters. 

at least as large as its width, the instantaneous ash residuals will contain information 
on the magnitude and position of the jet. These, combined with an initial guess of the 
mean strength, define the Gaussian velocity jet completely. A least squares fit to the 
Gaussian jet was performed on the a sh residuals for each subtrack and the parameters 
of the fit modified until a convergence criterion was satisfied. This method was applied 
to all tracks independently and yielded a mean ash field (Figure 4) that was consistent 
in position, but with a narrower and stronger GS than that obtained from hydrographic 
data. 

2*3 Combining the hydrographic and the synthetic means 

The Kelly /Gille/Qiu synthetic method only models the jet and cannot be used to recover 
the mean far from the GS itself. In other words, there is no information about the 
recirculation gyres that can be seen in the HydroBaae mean. To incorporate information 
away from the GS while retaining the narrower width of the jet, the simple Gaussian jet 
model was modified. If we assume that the difference in height between the synthetic and 
the HydroBase mean is due to the presence of northern and southern recirculation gyres 
in the latter, then the recirculation gyres can be modelled as a wide, slow Gaussian jet, 
whose height equals this difference (Figure 5). For this model, first mean height profiles 
of hydrographic data were obtained by splitting gridded HydroBaae data along satellite 
subtracki. An error function fit was performed on the profiles which yielded the function, 


8 



2 RECONSTRUCTING THE MEAN SEA SURFACE HEIGHT 



Figure 4: Synthetic aah mean obtained from TOPEX/Poseidon data. 

erf (a 1 , 03 , 03 ). where a lt as, aj are the amplitude, poiition and width parameter* of the 
fit. The recirculation wae then let to an error function er/(oj - a[, oj, 2 * 13 ), where aj i* 
the amplitude parameter of the lynthetic mean height profile. This recirculation error 
function was subtracted from the synthetic mean profiles to yield modified synthetic 
profiles. 

Finally, the modified synthetic profiles were gridded to 1° x I s and combined with the 
HydroBase mean using a spatially varying function. The profiles were first interpolated 
between tracks to retain the narrow jet structure of the Gulf Stream and then splined 
onto the I s x 1° grid using a bihamonic spline. The interpolation could be carried out 
since the mean height did not change much from track to track (except in the bifurcation 
region around 50°W). Figure 6 shows the spacing of the original tracks as thick solid 
lines, the interpolated ones are thin lines, and the open circles are the 1° x 1° grid 
point*. 

The weighting function used for combining the HydroBase mean with the gridded 
mean was calculated by low-pat* filtering the magnitude of the gradient of a$h of the 
•ynthetic mean (Figure 7). This procedure weighted the synthetic mean more heavily in 
regions of large velocities where the Gaussian model is expected to be more accurate, and 
the hydrographic mean more heavily away from the Gulf Stream, where the Gaussian 
model is not accurate. East of 40°W, where there was no synthetic mean data, the 
HydroBase mean wat used exclusively, whereai inside the 3000 meter Isobath, where no 
HydroBase data are available, the synthetic mean Was used. The resulting map is shown 
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2 RECONSTRUCTING THE MEAN SEA SURFACE HEIGHT 
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Figure 5: A comparison of total height (dashed line), HydroBase height (light line) and 
net height (heavy line) along eubtrack 228. 
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Figure 6: Heavy eubtrack lines mark the initial sampling pattern of the mean height, 
light lines show the interpolated eubtracks and circles indicate the final grid onto which 
it was mapped. 
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2 RECONSTRUCTING THE MEAN SEA SURFACE HEIGHT 



Figure 7: Weights used on the synthetic mean. 

in Figure 8. The combined mean field shows increased resolution over the hydrographic 
mean, and incorporates recirculation information which could not previously be obtained 
from the synthetic method. 

2.4 Obtaining Monthly Sea Surface Height Maps 

Mean ash field reconstructed in the previous step can now be added back to the tsh 
residuals to obtain monthly maps of the Gulf Stream and the North Atlantic Current. 
For this, the residuals need to be regridded in borh time and space. Assuming that the 
two are independent, the two-dimensional regridding can be broken up into two separate 
one-dimensional problems. The residuals are first averaged in time using a boxcar filter 
with a width of sixty days. This puts all cycles on a monthly spacing. The monthly 
cycles are then splined on to a i° by 1° grid using a two-dimensional biharaonic spline. 
They are then added to the means. See the Appendix for contour plots of the monthly 
maps so obtained. 

Monthly maps derived from the earlier Gtosat mission were also included in this 
study and covered a period from November 1988 to April 1989. Geosat Bubtr&ck data 
have an across-track spacing approximately half that of the TOFEX/Fostidon data, with 
a longer sampling period of 17 days. It was processed similarly but had an orbit cor- 
rection applied to it [2\ Because of smaller across-track spacing, it was gridded on to a 
0.5° x 0.5 s grid. Mean sea surface height was constructed by combining the synthetic 
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3 CALCULATING GULP STREAM INDICES 
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Figure 9: Error function fit to height along 50°W for November 1994. 
mean with historical hydrographic data [fj\, 

3 Calculating Gulf Stream Indices 

To characterize the annual cycle of meandering, as well as the strength and position of 
the Gulf Stream, a series of zonal indices were developed from monthly maps of *$h. 
These indices covered a region from 73° W to 50° W which was further broken down into 
two subregions: upstream and downstream of 63°W, the approximate location of the 
New England Seamounts, The calculation of each index is enumerated below: 


1. Transport: Height profiles along each degree of longitude are interpolated onto a 
0.125 degree grid, and error functions axe then fitted to them to estimate height 
difference across the jet. Figure 9 shows a fit for one profile. These height dif- 
ferences axe averaged over all longitudes to get the height difference index for the 
month. 

2. Position: If the path of the Gulf Stream it described by the zero contour of the 
height field, the mean position, P, is just a zonal average of the latitude of this 
zero contour. To get a robust estimate of path, the -0.1, 0 and 0.1 meter contours 
are averaged (Figure 10) before calculating P. 
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S CALCULATING GULF STREAM INDICES 



Figure 10: Path of the Gulf Stream jet (heavy black line) for November 1994, overlaid 
on ssh contours for that month. Slight differences in the path from the height contours 
are due to the different plotting packages used to create them. 

3. Curvature: Curvature of a function y = f(x) is given by: 

y* 

c = ITTTyTF (1) 

where the primes denote differentiation with respect to a. In this case, / is the 
path of the Gulf Stream as described in the calculation of P, z it the longitude of 
the path while y is the latitude. Curvature is calculated at every point along the 
path, and is then averaged across the path to yield the curvature index, C. 

4. Eastward Velocity: Eastward velocity is calculated along each degree of longitude. 
Since peak velocity estimates may be noisy, the velocity at a fixed width of one 
degree is used as the eastward velocity of the jet. This is averaged for all longitudes 
to yield the velocity index, V, 

The above four indices were calculated for both the TOPEX/Poseidoh at well as 
the G«o$at monthly maps. Gtosat derived monthly maps were decimated to a 1° x 
l e grid and GS indices were calculated similarly except for the position index, which 
was calulated using the center of the error function used for the ash difference. This 


14 



3 CALCULATING GULF STREAM INDICES 


was became of differences in the way the mean field was reconstructed in the two data 
sets as well as differences is spatial resolution. Figures 11 and 12 show the indices for 
TOPEX/Pdseidon and Gtotai for the upstream region (73°W to 64°W) , while Figures 13 
and 14 show the indices for the downstream region (63°W to 50 a W). Although there is no 
clear annual cycle apparent in the downstream indices, there is one >n the upstream ones. 
Putting together both time series and fitting an annual signal to them shows a maximum 
in position (corresponding to a minimum in curvature) in September with a minimum 
(and a corresponding maximum in curvature) in March. Transport is maximum in 
October and minimum in April. The amplitudes of the position and transport indices 
are 0.21° and 0.07m respectively. This implies that the Gulf Stream is stronger and 
follows a relatively straight and northeily path in fall, and then shifts to a southerly 
position with weaker transport accompanied by more meandering in spring. A detailed 
analysis of the index cycle may be found in [7]. 
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3 CALCULATING GULF STREAM INDICES 


TOPEX/Poseidon Upstream Indices 


a) Position 






Figure 11: Zonal indices for 73 a W to 64° W from TOPEX/Poitidon 
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Geosat Upstream Indices 




°'4dJ FTvttMJ JASCNDJFMAMJ JASONDJ FMA 
1986 1987 1988 1989 


d) Eastward Velocity 



Figure 12: Zonal indices for 73 9 W to 64®YV from Geosat. 
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S CALCULATING GULF STREAM INDICES 


a) Position 


TOPEX/Poseidon Downstream Indices 
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b) Transport 



d) Eastward Velocity 
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Figure 13: Zonal indices for 63°W to S0 S W from WPEXfPottldan. 
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APPENDIX 


Appendix 

Monthly maps of total stih derived from TtfPEX/Pastlddn altimeter data are included 
for the period from November 1992 to November 1994. The 3000 meter itobath (white 
contour line) it included in the maps for reference. 


21 








APPENDIX 


December 92 



80*W 70‘W 60‘W SO'W 40’W 


wnwMmu 


*0.8 -0.7 -0.6 -0.8 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

88H,m 


23 


APPENDIX 


January 93 


SO'N 


40'N 


30*N 



80'W 70' W 60'W 60*W 40*W 




■0.8 -0.7 *0.8 -0.6 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

SSH, m 


24 




APPENDIX 


February 93 



6Q‘W 70*W 6C‘W 50*W 40*W 




-O.fl -0.7 -0.8 -0.5 -0.4 -0,3 *0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.6 

88H,m 


25 


March 93 



26 




APPENDIX 


April 93 



80’W 70*W 60'W 50‘W 40‘W 


4 .1 sii'-PS, luVllMR 


•0.0 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

S8H, m 


27 


APPENDIX 


May 93 



80*W 70‘W 60 'W 60 4 W 40‘W 



-0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.6 0.6 0.7 0.8 


SSH.m 


28 




APPENDIX 


June 93 



80*W 70‘W 60‘W SOW 40*W 


J J&C&fclbNKMi 


-0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

SSH.m 


29 


APPENDIX 





APPENDIX 


August 93 



80*W 70*W 60'W 50*W 40’W 





-O.a -0.7 *0.6 -0.6 -0.4 -0,3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.6 0.6 0,7 0,8 

SSH.m 


31 


I 


APPENDIX 


September 93 



80*W 70'W 60'W 60‘W 40‘W 


' 1 


•0.8 -0.7 -0.6 *0.5 -0.4 -0.3 *0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

88H, m 


32 








APPENDIX 


October 93 



8Q*W 70'W 80*W SQ‘W 4Q‘W 


I !W! «BUffWa)ltefb.| 


-0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0,1 0.0 0.1 0.2 0.3 0.4 0.5 0.5 0.7 0.8 

S3H, m 


33 





APPENDIX 


November 93 



80'W 70*W 60‘W 5Q’W 40‘W 


rfSBMP 


•0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

SSH, m 


34 



APPENDIX 


December 93 



80*W 70‘W 60‘W {5Q'W 40‘W 



0.8 *0,7 *0.6 *0.6 *0.4 -0.3 *0.2 *0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


88H, m 


35 



APPENDIX 


January 94 



flO*W 70*W 60*W S0 4 W 40*W 


euvoFiw 


■0.8 -0.7 -0.6 -0.6 -0.4 -0.3 *0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

8SH, m 


36 







APPENDIX 


February 94 



80*W 70‘W 80‘W 50 "W 40*W 




i i i » i r i r i i i i 1 1 1 1 1 

-0.8 -07-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.6 

S8H, m 


37 





APPENDIX 


March 94 



60‘W 70*W 60'W 60*W 40‘W 


im ikm# «n9»Kasanr 


•0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

SSH, m 


38 



APPENDIX 


April 94 



60*W 70'W 60*W 50*W 40‘W 





■0.8 -0.7 -0.6 -0.5 -0.4 *0.3 -0,2 -0.1 0.0 0.1 0.2 0.3 0.4 0.6 0.6 0.7 0.6 

MH, m 


39 



APPENDIX 


May 94 



- 0.6 - 0.7 * 0.8 - 0.5 * 0.4 - 0.3 - 0.2 - 0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


«SH,m 


40 



APPENDIX 


June 94 



80'W 70‘W 60’W 50 ’W 40*W 


! v&nRaiNMM 


-0.8 -0.7 -0.0 -0.6 -0.4 -0,3 *0,2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.8 0.7 0.8 

SSH.m 


41 



APPENDIX 


July 94 



80'W 70*W 60’W 50*W 40‘W 




I I « T I 1 T I I I I | | “I 1 1 

•0.8 -0.7 -0.6 -0.5 *0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.6 0.6 0.7 0.8 

SSH.m 


42 




APPENDIX 


August 94 



80*W 70*W 60‘W 50*W 40*W 


p" m m 


i * i i i i 1 i i i 1 1 1 1 r 1 1 

■0.8 -0.7 *0.0 *0.5 -0.4 -0.3 -0.2 *0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

8SH, m 


43 



APPENDIX 


September 94 



• 0.8 ” 0.7 - 0.8 - 0.5 - 0.4 - 0.3 - 0,2 - 0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


SSH.m 


44 




APPENDIX 


October 94 



80*W 70'W 0Q*W 50’W 40‘W 



■0.8 *0.7 *0.0 *0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


S8H, m 


46 


APPENDIX 


November 94 



80'W 70‘W 80'W 50'W 40‘W 


* f f M* MMtMBHDU 


■0.8 -0.7 -0.6 *0.5 -0.4 -0.3 -0.2 -0,1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

SSH, m 


46 





DOCUMENT LIBRARY 

Dirtribution Lift for Technical Report Exchange * February 199b 


University of California, San Diego 
SIO Library 0175C 
9500 Gilman Drive 
La Jolla, CA 92093-0175 

Hancock Library of Biology Sc Oceanography 

Alan Hancock Laboratory 

University of Southern California 

University Park 

Los Angeles, CA 90089-0371 

Gifts Sc Exchanges 
library 

Bedford Institute of Oceanography 
P.O. Box 1006 

Dartmouth, NS, B2Y 4A2, CANADA 

Commander 
Intematl~4ial Ice Patrol 
1082 Shennecossett Road 
Groton, CT 06340-6095 

NOAA/EDIS Miami Library Center 
4301 Rickenbacker Causeway 
Miami, FL 33149 

Research Library 
U.S. Army Corps of Engineers 
Waterways Experiment Station 
3909 Halls Ferry Road 
Vicksburg, MS 39180-6199 

Institute of Geophysics 
University of Hawaii 
Library Room 252 
2525 Correa Road 
Honolulu, HI 96822 

Marine Resources Information Center 

Building E38-320 

MIT 

Cambridge, MA 02139 
Library 

Lamont-Doherty Geol jical Observatory 
Columbia University 
Palisades, NY zl0964 

Library 

Serials Department 
Oregon State University 
Corvallis, OR 97331 

Pell Marine Science Library 
University of Rhode Island 
Narragansett Bay Campus 
Narragansett, R1 02882 


Working Collection 
Texas A&M University 
Dept of Oceanography 
College Station, TX 77843 

Flsheries-Oceanography Library 
151 Oceanography Teaching Bldg. 
University of Washington 
Seattle, WA 98195 

Library 

R.S.M.A.S. 

University of Miami 
4600 Rickenbacker Causeway 
Miami, FL 33149 

Maury Oceanographic Library 
Naval Oceanographic Office 
Building 1003 South 
1002 Batch Blvd. 

Stennis Space Center, MS, 39522-5001 
Library 

Institute of Ocean Sciences 
P.O. Box 6000 
Sidney, B.C.V8L4B2 
CANADA 

National Oceanographic Library 
Southampton Oceanography Centre 
European Way 
Southampton SO 14 3ZH 
UK 

The Librarian 

CSIRO Marine Laboratories 
G.P.O. Box 1536 
Hobart, Tasmania 
AUSTRALIA 7001 

Library 

Proudmah Oceanographic Laboratory 

Bldston Observatory 

Birkenhead 

Merseyside L43 7 RA 

UNITED KINGDOM 

IFREMER 
Centre de Brest 

Service Documentation • Publications 
BP 70 29280 PLOUZANE 
FRANCE 


IMPORT DOCUMnff ATWN 
PAQK 


WHOI-97-06 


$ of Sea Surface Height in the North Atlantic and Zonal Indices I &RWC Juoe 1997 
for the Gulf Stream Using TOPEX/Poseidon Altimeter Data 


for the Gulf Stream Using TOPEX/Poseidon 


Sandipa Singh and Kathryn A. Kelly 


OhMHMkMNMIlMdAdiVM 

Wf nWflf VW 

Woods Hole Oceanographic institution 
Woods Hole, Massachusetts 02543 



ii.OeittMeo erOnotfOi No. 
(6) NAOW-l<566 
NAGW-4806 


1litp0nHnn|Or|nMiMNffWin0A<Miii 

NASA 


IS. Type et mtMfl a PwtM Oovme 
Technical Report 


This report should be died as: Woods Hole OceanOg. Inst Tech. Rept, WHOI-97-06. 


ia AMeaet ammc see wares) 

Monthly Mips of sea surface height are constructed for die North Atlantic Ocean using TOPEX/Ponidcn altimeter data. 
Mean sea surface height Is reconstructed using a weighted combination of historical, hydrographic data and a synthetic 
mean obtained by fitting a Gaussian model of the Gulf Stream jet to altimeter data. The resultant mean shows increased 
resolution over the hydrographic mean, and incorporates recirculation information that is absent in the synthetic mean. 
Monthly maps, obtained by adding the mean field to aldmeier sea surface height residuals, are used to derive a set Of zonal 
indices that describe the annual cycle of meandering as well as position and strength of the Gulf Stream. 


17. Mwnthi Annual* aMaertetore 
altimeter 

Gulf Stream 
monthly maps 


A MemrteratOptivrttote Tims 


AOOOATtnewarMp 


i| A Is. 

W. MfMHR flf MWMni 


Approved for public rclcsic; distribution unlimited 


(St* ANttZtS. 1 l) 



«Ml^rud«0nign%ywN 


OPTIONAL POfIM ITS ( 4 * 77 ) 

(Form*<ty MTIS.J0) 
Otps/inwrt d Cerrmwo* 








